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ABSTRACT 



High-altitude pulmonary hypertension (HAPH) is a consequence of chronic alveolar hypoxia, leading to hypoxic vasoconstriction 
and remodeling of the pulmonary circulation. Brisket disease in cattle is a naturally occurring animal model of hypoxic pulmonary 
hypertension. Genetically susceptible cattle develop severe pulmonary hypertension and right heart failure at altitudes >7,000 ft. 
No information currently exists regarding the identity of the pathways and gene(s) responsible for HAPH or influencing severity. 
We hypothesized that initial insights into the pathogenesis of the disease could be discovered by a strategy of (1 ) sequencing of 
functional candidates revealed by single nucleotide polymorphism (SNP) analysis and (2) gene expression profiling of affected 
cattle compared with altitude-matched normal controls, with gene set enrichment analysis (GSEA) and Ingenuity pathway analysis 
(IPA). We isolated blood from a single herd of Black Angus cattle of both genders, aged 12-18 months, by jugular vein puncture. 
Mean pulmonary arterial pressures were 85.6±13 mmHg STD in the 10 affected and 35.3±1.2 mmHg STD in the 10 resistant cattle, 
P<0.001. From peripheral blood mononuclear cells, DNAwas hybridized to an Affymetrix 10K Gene Chip SNP, and RNAwas used 
to probe an Affymetrix Bovine genome array. SNP loci were remapped using the Btau 4.0 bovine genome assembly. mRNAdata 
was analyzed by the Partek software package to identify sets of genes with an expression that was statistically different between 
the two groups. GSEA and I PA were conducted on the refined expression data to identify key cellular pathways and to generate 
networks and conduct functional analyses of the pathways and networks. Ten SNPs were identified by allelelic association and 
four candidate genes were sequenced in the cohort. Neither endothelial nitric oxide synthetase, NADH dehydrogenase, TG- 
interacting factor-2 nor BMPR2 were different among affected and resistant cattle. A 60-gene mRNA signature was identified that 
differentiated affected from unaffected cattle. Forty-six genes were overexpressed in the affected and 14 genes were downregulated 
in the affected cattle by at least 20%. GSEA and Ingenuity analysis identified respiratory diseases, inflammatory diseases and 
pathways as the top diseases and disorders (P<5.14*10" 14 ), cell development and cell signaling as the top cellular functions 
(P<1.20x10 08 ), and IL6, TREM, PPAR, NFkB cell signaling (P<8.69*10-° 9 ) as the top canonical pathways associated with this 
gene signature. This study provides insights into differences in RNA expression in HAPH at a molecular level, and eliminates 
four functional gene candidates. Further studies are needed to validate and refine these preliminary findings and to determine 
the role of transcribed genes in the development of HAPH. 
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INTRODUCTION 

High mountain disease (brisket disease] is right heart 
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failure due to hypoxic pulmonary hypertension in cattle 
residing at high altitudes. 1131 Hypoxia is the most potent 
stimulus for pulmonary hypertension, and the hypoxia of 
high altitude (>7,000 ft] is a well known cause. 14 81 Some 
cattle [Bos taurus) possess a heritable susceptibility to 
severe high-altitude pulmonary hypertension (HAPH]. [912] 
While most cattle thrive at high altitudes, susceptible 
cattle develop pulmonary hypertension that is sufficient 
to cause right heart failure, edema of the brisket and 
death. Experiments conducted by Grover and Reeves 
in cattle susceptible and resistant to HAPH suggest an 
autosomal-dominant mode of inheritance, or inheritance 
related to a few major genes. [1115] No information 
currently exists on the identity of these genes. HAPH in 
cattle occurs in about 20% [10-50%] of animals brought 
to high altitudes (>7,000 ft] to replenish herds, [81416] 
suggesting that the gene of interest may be a relatively 
common polymorphism only expressed with the stimulus 
of altitude hypoxia. There is no known phenotype of the 
gene at low altitude. 

We hypothesized that the insights into pathogenesis of 
brisket disease may be gained through two approaches: 
(1] allelic association analysis of single nucleotide 
polymorphism (SNP] differences; and [2] a comparative 
gene expression analysis of hypertensive versus resistant 
cattle. Therefore, we analyzed DNA and RNA from 20 cattle 
of the same herd; 10 with severe pulmonary hypertension 
and 10 with normal pressures. We used an Affymetrix 
10K Bovine Gene SNP array and sequenced four genes of 
interest. Expression data were then analyzed for genes 
and canonical pathways that were statistically different 
between the two groups. Gene set enrichment and 
Ingenuity map analyses were performed to further refine 
the data. 116191 Collectively, the results of this study provide 
molecular and cellular observations in possible genes of 
interest in HAPH leading to brisket disease. 

MATERIALS AND METHODS 

Selection of study animals and blood samples 

We obtained blood from a single herd of black Angus 
cattle of both genders, aged 12-18 months, residing at 
8,500 feet, by jugular vein puncture on the same day of 
right heart catheterization screening for HAPH by one 
of us (Timothy N. Holt, hereinafter TH]. The project was 
approved by the Vanderbilt Medical Center IACUC. Right 
heart catheterization is standard practice at altitude to 
discover hypertensive animals so that they can be shipped 
down to a low altitude. 

DNA and RNA preparation and analysis 

DNA was isolated from whole blood using QIAmp DNA 
mini-kit as directed by the manufacturer's instructions 



(Qiagen, Valencia, CA, USA] and subsequently quantified 
using a spectrophotometer. Total cellular RNA was 
prepared from blood on silica-gel membranes with on- 
column DNAse treatment (RNeasy Total RNA Isolation 
Kit; Qiagen]. RNA integrity was assessed via agarose 
gel electrophoresis and samples showing evidence of 
degradation were discarded. 

Affymetrix bovine 10K SNP genotyping and 
bioinformatic analysis 

For SNP comparison, genomic DNA was extracted as 
described above and hybridized to Affymetrix - GeneChip 
Bovine Mapping 10 K SNP chip according to manufacturer's 
instructions (Affymetrix, Santa Clara, Calif., USA]. 
Following bovine genotype collection, we remapped 
all of the SNP loci represented on the bovine 10K array 
using the Batch query within dbSNP online (NCBI] in 
conjunction with the latest bovine genome assembly 
(Btau 4.0]. In addition, SNP loci present on the array 
more than once (duplications] were identified, and 
all duplications were removed. Likewise, SNPs with 
multiple chromosomal assignments were also flagged 
and subsequently removed. At the conclusion of our 
bioinformatic analyses, 8,011 unique bovine SNPs with 
single chromosomal assignments and updated bovine 
genome annotation were available for genome-wide 
association (GWA] analysis. 

Statistical approach to SNP microarray 

Because the precise mode of inheritance related to 
HAPH in cattle has not been fully elucidated, we used 
several standard case-control GWA approaches in a 
preliminary effort to identify new potential candidate 
genes. All statistical tests were performed using formulae 
incorporated within the software program HelixTree 
6.4.2 (Golden Helix Inc., Bozeman, MT, USA]. In the first 
approach, we performed a simple allelic association test 
across all unique SNP loci for 10 cattle exhibiting severe 
HAPH (n=20 alleles] and 10 cattle defined as "resistant" 
at altitude (n=20 alleles]. Thereafter, we also performed 
standard genotypic association tests in an effort to 
evaluate whether SNPs implicated by a simple allelic test 
would systematically assemble into disparate genotypic 
classes for cattle with severe HAPH, as compared with 
altitude-matched HAPH-resistant controls. Raw P-values 
for all tests were corrected by either full-scan permutation 
and/or the Bonferroni procedure implemented within 
HelixTree 6.4.2. 

Sequencing SNP candidates 

TG-interacting factor-2 (TGIF2], endothelial nitric 
oxide synthetase (eNOS] and NADH coding regions 
were sequenced using a cDNA-based approach. Briefly, 
we performed first strand cDNA synthesis using 
the Superscript First-Strand System (Invitrogen Life 
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Technologies, Carlsbad, CA, USA] with 3 mg of total RNA 
and an oligo-dT primer. We then used one-tenth volume 
of the first strand reaction as a template for polymerase 
chain reaction (PCR] amplification using the Elongase 
Amplification System (Invitrogen Life Technologies]. The 
forward and reverse primers corresponded to sequences 
located in the 5' and 3' untranslated mRNA sequence 
of BMPR2, eNOS and TGIF-2 (primer sequences and 
conditions available upon request]. We then carried out 
individual PCR reactions and visualized all amplicons 
on a 2% agarose gel. Thereafter, we used 2 ul ExoSAP-it 
[USB] per 5 ul of PCR reaction (37°C for 15 min, followed 
by 80°C for 15 min] to eliminate residual primer and 
dNTPs. Sequencing was accomplished using the Big 
Dye Terminator cycle sequencing technology (Applied 
Biosystems, Carlsbad, Calif., USA] in 10 ul reactions 
containing 2 ul Big Dye, 2 ul Big Dye Buffer, 50 ng template 
and 2.5 ul of primer. Thermal cycling parameters were as 
follows: 95°C for 5 min; 30 cycles at 95°C for 30 s, 55°C 
for 10 s and 60°C for 4 min. After ethanol precipitation, 
the products were resuspended in 10 ul formamide 
and analyzed on a 3100 Genetic Analyzer (Applied 
Biosystems]. 

Microarray and gene set enrichment analysis 
analyses 

Pulmonary tissue was not available from these animals. 
It is also possible that the lung tissue expression 
signature might be overwhelmed by adaptive cellular 
responses. This would then make the interpretation of 
the expression data particularly susceptible to false- 
positive and false-negatives. 116191 Our solution to this was 
to avoid disease effector cells and to use peripheral blood 
mononuclear cells (PBMCs] in our studies. Peripheral 
blood mononuclear cells are genetically the same as lung 
tissue intraindivi duals, and are genetically different among 
hypertensive and resistant cattle and are likely to be free 
of end-stage disease effects. RNA was extracted from 
PBMCs. After the RNA was isolated, it was used to probe 
the Affymetrix GeneChip Bovine Genome Array (900561]. 
Raw data were then analyzed by the Partek software 
package at the Vanderbilt Functional Genomics Shared 
Resource to produce expression intensities for each probe 
set. To study pathway-level differences between groups, 
GSEA was conducted on the microarray gene expression 
data. 111 We used the Gene Ontology Biological Process 
(c5; from the MSigDB database] as our gene set database. 

In silico functional analysis HPAH signature 
using ingenuity pathways analysis 

We used IPA in the Ingenuity System to generate networks 
and conduct functional analyses of the HAPH signature. 
A data set containing gene identifiers was uploaded into 
the application. These genes, called focus genes, were 



used to query a global molecular network developed 
from information contained in the Ingenuity Pathways 
Knowledge Base. IPA generates models of gene interactions 
called networks that are presented graphically to show 
relationships between genes and the pathways they 
regulate. These networks are ranked according to a score 
calculated via a right-tailed Fisher's exact test. In network 
graph, proteins encoded by genes are represented as 
nodes and their relationships as edges (links]. All edges 
are supported by references from the literature. The 
functional analysis of a network identified the biological 
functions and/or diseases that were most significant to 
the genes in the network. The network genes associated 
with biological functions and/ or diseases in the Ingenuity 
Pathways Knowledge Base were considered for further 
analysis. Fisher's exact test was used to calculate a P-value, 
determining the probability that each biological function 
and/or disease assigned to a network is by chance alone. 
Canonical pathways analysis identified the pathways 
from the IPA library of canonical pathways that were 
most significant to the data set. The significance of the 
association between the data set and a canonical pathway 
was measured in two ways. 111 A ratio of the number of genes 
from the data set that map to the pathway divided by the 
total number of genes that map to the canonical pathway 
is displayed. 121 Fischer's exact test was used to calculate a 
P-value determining the probability that the association 
between the genes in the dataset and the canonical 
pathway is explained by chance alone. 

RESULTS 

Study animals 

Pulmonary arterial (PA] pressures were measured 
by jugular vein puncture at the time of right heart 
catheterization. Mean PA pressures were 86.10±5 mmHg 
STD in the affected and 31.20±0.7 mmHg STD in the 
resistant cattle, P<0.0001, by two-tailed Upaired t test 
(Fig. 1]. A mean PA pressure of 30-40 mmHg is normal for 
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Figure 1 : Mean pulmonary arterial pressures in the 1 0 cattle with pulmonary 
hypertension at altitude versus 1 0 unaffected at the same altitude. 
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cattle residing at altitude for at least 12 months. 13 ' 81 A mean 
PA pressure greater than 49 mmHg denotes a high risk for 
development of brisket disease. 181 Inclusion criteria were 
animals with a mean pulmonary artery pressure (mPAP] 
of 72-116 mmHg, thereby considered "affected," and 
those with a mPAP of 32-37 mmHg, thereby considered 
"resistant" at altitude. 

GWA using the affymetrix geneChip bovine 
mapping 1 0K array 

Among the bovine cohort investigated, 6,344 of the 
8,011 unique bovine SNPs with single chromosomal 
assignments and updated bovine genome annotation 
were diallelic, thereby making them potentially 
informative for subsequent analyses. The most robust 
allelic association was detected on BTA10 (0.0000028 
< p raw ^ 0.0000451], with weaker signals detected on 
additional bovine chromosomes. As expected, full-scan 
permutation and/ or Bonferroni correction for multiple 
testing revealed a much weaker BTA10 signal for both 
approximate and exact tests (0.0476 < p corrected ^ 0.2856], 
respectively. Summary data for the top 10 bovine SNPs 
implicated by an allelic GWA test are presented in Table 1. 
Additionally, bovine genes proximal to each SNP are also 
given, with relevant human-bovine comparative data. 

Examination of relevant genotypes determined for the 
top 10 SNP loci elucidated by an allelic association test 
revealed six corresponding loci whereby severe bovine 
HAPH cases were fixed for a genotype underrepresented 
in our altitude-matched HAPH-resistant controls (Table 1]. 
However, genotypic association tests for all variable loci 
(n=6344] revealed no statistically significant associations 
after correction for multiple testing. Relevant bovine 
genotypic data and corresponding distributions are 
summarized in (Table 1]. 



Because of its role in human PAH [20 " 231 and its place in the 
TGF-b family of growth and repair receptors, the bovine 
BMPR2 gene was an obvious initial functional candidate for 
HAPH. We screened 10 severely affected and 10 unaffected 
cattle with a set of bovine microsatellites proximal to 
bovine BMPR2 to evaluate whether local genetic variation 
exhibited any evidence of segregation with the brisket 
disease phenotype (HAPH]. This microsatellite analysis 
revealed a total of seven BMPR2 alleles, with no significant 
differences in the allelic distributions between the HAPH- 
affected and -resistant cattle, as determined by Fisher's 
exact test. No microsatellite allele was found uniformly 
among severely hypertensive cattle or resistant cattle thus 
providing no evidence of linkage disequilibrium between 
any particular allele and the HAPH phenotype. These data 
suggest, but do not prove, that variation within bovine 
BMPR2 is not responsible for modulating HAPH in cattle. 

Two other high-priority functional candidates were 
screened by cDNA-based sequence analysis, including eNOS 
and TGIF2. [24 ' 25] Low NO production has been implicated 
in the pathogenesis of high-altitude pulmonary edema 
(HAPE], [26 271 while TGIF2 is a homeobox transcriptional 
repressor and an important downregulator of BMPR2 and 
TGF-b signaling. 127 281 Sequence comparison of the eNOS 
and TGIF2 coding regions of HAPH-affected and -resistant 
cattle did not result in the identification of any statistically 
significant variants predictive of HAPH in our samples. 

Statistical analysis of microarray data reveals 
altered expression in HAPH cattle 

Gene expression profiles were determined by 
oligonucleotide microarray analyses using total RNAs 
from the PBMCs cells as described in the Methods section. 
In order to identify candidate genes involved in HAPH or 
protection from HAPH, we used a hierarchical clustering 
analysis to analyze the expression signature from the 



Table 1 : Top six bovine SNPs implicated by allelic association tests for a small bovine cohort consisting of 
severe HAPH cases and altitude-matched HAPH-resistant controls 



Bovine 


Bovine 


Proximal bovine 


Human 


Human Mb 


Relevant 


SNP a 


chrom. b 


gene(s) c 


chrom. d 


range e 


human gene' 


rs29024744 


BTA10 


NRXNIII/DI02 


14 


79,733,622- 


Deiodinase 


rs29025872 


BTA24 


BRUNOL4 


18 


79,748,2 
33,077,828- 


RNA binding 


rs29024519 


BTA29 


FOXR2 


X 


33,391,941 
55,666,558- 


protein 
Forkhead 


rs29016420 


BTA1 


MYH15 


3 


55,668,483 
109,581,906- 


box R2 
Myosin heavy 


rs29024944 


BTA24 


NDUFV2/LOC781824 


18 


109,730,859 
9,092,725- 


chain 15 
NADH 


rs29013210 


BTA3 


LOC100138920/FKBP1A 


20 


9,124,313 
1,297,622- 
1,321,753 


dehydrogenase 
FK binding 
protein 1A 



SNP: single nucleotide polymorphism; HAPH: high-altitude pulmonary hypertension; NADH: Nicotinamide adenine dinucleotide phosphate 
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PBMCs. A combination of relative level (fold-change) 
of expression and statistical significance (Student's 
t-test] was used to distinguish these genes. This analysis 
identified 675 genes that were significantly upregulated 
and 320 genes that were downregulated (the HAPH 
signature] in affected animals. The top 15 candidates are 
listed in Table 2. Interestingly the gene with the highest 
expression in affected animals is important for antigen 
presentation while the gene with the lowest expression 
is important for vessel elasticity as well as regulation of 
the TGF-beta pathway (Table 2). Genes related to disease 
process are shown in (Table 3). 

DISCUSSION 

HAPH with right heart failure was first described by Glover 
and Newsome in 1915 as "Brisket disease: Dropsy of high 
altitude." 111 The link of altitude to hypoxic pulmonary 
hypertension was made after the discovery of hypoxic 
pulmonary vasoconstriction by von Euler and Liljestrand 
in 1946. [28] The exact mechanism causing acute hypoxic 
vasoconstriction remains elusive, although much is 



known about modifying influences. [29 31] Chronic hypoxia 
is a powerful stimulus for pulmonary hypertension, which 
involves not only vasoconstriction but also remodeling 1321 
of the pulmonary arteries. K-channel function has a 
central role, as does oxidant stress, 133341 with Rho-Rho 
kinases involved in transduction of the hypoxic pressor 



response 



[35-39] 



To date, no major genetic variants have been identified that 
modify the strength of the hypoxic pressor response. [37] 
Some observations have been made regarding genetic 
associations with pulmonary responses to altitude in 
humans. 130 401 For example, eNOS and tyrosine kinase 
gene variants were overrepresented in a Japanese cohort 
of individuals susceptible to HAPE as compared with 
healthy controls at altitude. 127 401 Eddahibi and colleagues 
detected an association between the LL genotype of the 
human SLC6A4 (formerly SERT) with pulmonary artery 
pressure in a cohort of patients with chronic obstructive 
pulmonary disease as compared with controls. 1411 Recent 
studies implicate HIF2a polymorphism in Tibetans versus 
Han as a possible difference in successful high-altitude 



Table 2: The top 15 up and down genes expressed in peripheral blood mononuclear cells comparing high- 
altitude pulmonary hypertension with resistant cattle at 8,500 ft altitude 



Name 


Fold-change 


Molecular function 


Functions in cells 


HLA-DQA2 


+4.95 


MHC Class II receptor 


Antigen processing and presentation 


DNER 


+2.97 


Notch binding and calcium ion binding 


Differentiation, notch signaling 


SLC28A3 


+ 2.88 


Sodium symporter 


Vascular tone 


RYR1 


+ 2.43 


Voltage-gated calcium channel 


Response to hypoxia, calcium transport 


FCAR 


+ 2.36 


IgA binding 


Immune response, respiratory burst 


RAB3IP 


+ 2.25 


GTPase binding 


Golgi to plasma membrane transport 


ARRDC4 


+ 2.21 


Unknown binds MLX, MLXIP 


Unknown 


UPP1 


+ 2.17 


Uridine phosphorylase 


Growth 


GPD2 


+2.10 


Glycerol-3-phosphate dehydrogenase 


NADH metabolic process, oxidation 






[NAD+] activity 


reduction 


TGM3 


+2.05 


Calcium ion biding, GTPase activity 


Differentiation, protein modification 


EMR1 


+2.01 


G protein-coupled receptor 


Proliferation, adhesion, calcium response 


HIST1H1C 


+ 1.93 


DNA binding 


DNA binding 


PTX3 


+ 1.93 


Zymosan binding 


Inflammatory response, +NOS process 


NFIL3 


+ 1.78 


Transcription factor 


Immune function, particularly T cells 


DEFB4A 


1.77 


Chemotaxis, immune response 


Immune response 


FBN1 


-2.03 


Extracellular matrix, regulation of TGF-beta 


Muscle and heart development 


KIR3DL1 


-1.6 


Inhibitory MHC Class I receptor 


Regulation of immune response 


BOLA 


-1.2 


Class I MHC 


Antigen processing, immune response 


CRIP3 


-1.15 


Zinc and metal ion binding 


T cell function and proliferation 


CLGN 


-1.13 


Calcium ion binding, unfolded protein response 


Endoplasmic reticulum chaperone protein 


IKZF3 


-1.01 


Zinc-binding transcription factor 


Mesoderm development, lymphocyte 








development 


TCRB 


-0.97 


T-cell receptor 


Immune response 


RAB30 


-0.89 


GTP binding, GTPase activity 


GTP signaling, protein transport 


EBF1 


-0.82 


Transcription factor 


Apoptosis 


PHACTR1 


-0.78 


Actin binding 


Actin function 


CD8A 


-0.75 


MHC Class I 


Immunity, calcium-mediated signaling 


BCL11A 


-0.74 


Transcription factor 


Apoptosis, differentiation 


TCRA 


-0.67 


T-cell receptor 


Immune response 


MY05C 


-0.66 


Myosin complex 


Muscle function 


VEGFB 


-0.65 


Endothelial growth factor 


Angiogenesis, positive regulation of 








endothelial function 
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Table 3: Top 10 disease processes associated with high-altitude pulmonary hypertension in cattle 



Category 



P-value Genes 



Respiratory disease 



Inflammatory 
response 



Inflammatory 
disease 



Connective tissue 
disorders 



Skeletal and 
muscular disorders 



Immunological 
disease 



Genetic disorder 



Hematological 
disease 

Cardiovascular 
disease 



Metabolic disease 



1.61E-11- 
7.01E-03 

1.02E-10- 
7.81E-03 



1.66E-08- 
7.01E-03 



1.07E-07- 
4.28E-03 



1.07E-07- 
6.72E-03 



2.05E-07- 
5.92E-03 



1.33E-06- 
6.72E-03 



2.51E-06- 

6.05E-05- 
6.83E-03 



1.12E- 



HIST1H1C, SOCS3, TREM1, ACVRL1, ARG2, CD164, HRH1, HMOX1, NFIL3, GPR77, 
CXCR2, F5, MXD1, RYR1, FOSL2, LTBR, PDXK, FCGR3A, PTX3, IL15, IL1R1, TLR4, 
CSF2RB, HP, NCF1, S100A9, IL1RN, CEBPD, FBN1, IL1B, CD14, IL2RA, CHI3L1 
FCAR, KIR3DL1, SOCS3, TREM1, IGSF6, OLFM4, ETS2, CASP4, BCL6, IL1R2, CD164, 
HMOX1, HRH1, NFIL3, CD47, GPR77, SOD2, HLA-A, CXCR2, F5, S100A8, LTBR, CCL16, 
FCGR3A, TNFSF13B, PTX3, CCR1, MGST1, ATF3, IL15, ALOX5AP, IL1R1, CSF2RB, TLR4, 
AQP9, NCF1, HP, S100A9, DUSP1, IL1RN, CEBPD, IL1B, CD14, IL2RA, SIRPA 
SOCS3, KIR3DL1, DPYD, ARG2, MAPK13, BCL6, IL1R2, HRH1, GPR77, SOD2, GPD2, 
FOSL2, DYSF, LTBR, SLC04C1, TNFSF13B, ACSL6, MYOF, STX2, ACVR1B, CSF2RB, 
CLIC2, NCF1, AQP9, IL1RN, DUSP1, FBN1, CD14, RNF149, HLA-DQA2, SIRPA, MCTP2, 
TREM1, RBM47, CASP4, HMOX1, WDFY3, HLA-A, CXCR2, F5, WLS, MXD1, S100A8, 
C9ORF150, FCGR3A, CCR1, KCNK17, DCK, DNER, IL15, ALOX5AP, IL1R1, TLR4, HRSP12, 
P2RY13, HP, S100A9, FAM129A, IL1B, IL2RA, C150RF48, EMR1, CHI3L1, TCF7L2, MAOA 
SOCS3, DPYD, MAPK13, BCL6, IL1R2, HMOX1, SOD2, HLA-A, GPD2, CXCR2, ANTXR2, 
S100A8, LTBR, TNFSF13B, SLC04C1, FCGR3A, CCR1, KCNK17, DNER, IL15, ALOX5AP, 
IL1R1, MYOF, ACVR1B, TLR4, CSF2RB, P2RY13, NCF1, HP, CLIC2, AQP9, S100A9, 
DUSP1, IL1RN, FAM129A, IL1B, FBN1, IL2RA, RNF149, CHI3L1, HLA-DQA2, SIRPA, 
TCF7L2 

SOCS3, DPYD, SGK1, MAPK13, BCL6, NPL, IL1R2, HRH1, GPCPD1, SOD2, GPD2, PDXK, 
DYSF, LTBR, SLC04C1, TNFSF13B, MYOF, ACVR1B, KAT2B, CSF2RB, NCF1, AQP9, CLIC2, 
IL1RN, DUSP1, CD14, FBN1, RNF149, TBXAS1, MTDH, HLA-DQA2, SIRPA, MCTP2, ETS2, 
GNG7, HMOX1, HLA-A, CXCR2, MXD1, F5, SNX10, RYR1, S100A8, FCGR3A, CCR1, 
KCNK17, DCK, DNER, IL15, ALOX5AP, IL1R1, TLR4, P2RY13, HP, S100A9, FAM129A, 
IL1B, IL2RA, CHI3L1, TCF7L2, 

SOCS3, KIR3DL1, DPYD, MAPK13, FNDC3B, BCL6, IL1R2, HRH1, GPD2, LTBR, DYSF, 
SLC04C1, TNFSF13B, MYOF, ACVR1B, CSF2RB, NCF1, AQP9, DUSP1, IL1RN, CD14, 
FBN1, RNF149, HLA-DQA2, ELL2, SIRPA, TREM1, CASP4, ID1, HMOX1, CD47, HLA-A, 
WDFY3, CXCR2, S100A8, FCGR3A, CCR1, OAS1, KCNK17, DCK, DNER, PSTPIP2, IL15, 
IL1R1, TLR4, HRSP12, P2RY13, HP, S100A9, FAM129A, IL1B, IL2RA, CHI3L1, TCF7L2 
SOCS3, DPYD, KCNJ2, SGK1, KLF6, ACVRL1, ARG2, FNDC3B, BCL6, NPL, IL1R2, 
GPCPD1, HRH1, SOD2, GPD2, FOSL2, DYSF, ALPK1, PDXK, LTBR, TNFSF13B, SLC04C1, 
MGST1, ATF3, ACSL6, MYOF, VLDLR, STX2, SERINC2, KAT2B, CSF2RB, NCF1, IL1RN, 
DUSP1, FBN1, CD14, RFX2, MTDH, TBXAS1, HLA-DQA2, SIRPA, MCTP2, LITAF, TREM1, 
RBM47, ETS2, CASP4, VAT1L, EIF4E, GNG7, HMOX1, ID1, CD47, HLA-A, WDFY3, CXCR2, 
MXD1, F5, WLS, ADAM 19, SNX10, ANTXR2, RYR1, S100A8, C9ORF150, MARCKS, 
AGPAT9, RNF144B, RAB8B, FCGR3A, CCR1, MOCS1, KCNK17, DCK, DNER, IL15, DGAT2, 
ALOX5AP, IL1R1, PPP1R3B, MARCKSL1, TLR4, HP, S100A9, FAM129A, IL1B, IL2RA, 
EMR1, C150RF48, CHI3L1, TCF7L2, ISG20, MAOA 

SOCS3, TREM1, MAPK13, BCL6, HMOX1, HRH1, CD47, SOD2, GPD2, F5, RYR1, FOSL2, 
LTBR, FCGR3A, ATF3, DCK, IL15, TLR4, CSF2RB, S100A9, IL1RN, IL1B, CD14, IL2RA, SIRPA 
SOCS3, DPYD, KCNJ2, KLF6, ACVRL1, ARG2, BCL6, NPL, HRH1, SOD2, GPD2, ALPK1, 
DYSF, PTX3, MYOF, VLDLR, KLF4, CSF2RB, NCF1, AQP9, IL1RN, CD14, FBN1, SIRPA, 
RBM47, ETS2, VAT1L, HMOX1, CD47, F5, ANTXR2, ADAM19, RYR1, S100A8, KCNK17, 
PSTPIP2, IL15, DGAT2, ALOX5AP, IL1R1, TLR4, MARCKSL1, HP, S100A9, IL1B, IL2RA, 
CHI3L1, ISG20, MAOA 

ID1, HMOX1, CD47, SOD2, IL1RN, DUSP1, CEBPD, CD14, IL1B, BCL6, KLF4 



adaptation. [42] Many other signaling pathways possess 
potential candidate genes.' 431 

Relevant to the present study, the risk of acquiring brisket 
disease in low-altitude cattle brought to very high altitude 
(10,000 ft) is about 10-50%.P 144 i Notably, the prevalence of 



brisket disease among herds chronically residing and bred 
above 2100 m (c. 7,000 ft] has been reduced to about 1%, 
presumably because of selective loss of susceptible cattle 
to heart failure and to removal from herds after HAPH has 
been identified. 131 Previously, Reeves, Weir and Grover bred 
a cohort of cattle with HAPH whose mean PA pressure was 
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50mmHg (n=8) and a group of resistant animals whose 
mean PA was 29mmHg (n=ll), all residing at 10,000 ft. [1114] 
The 19 cattle were initially taken down to an altitude of 
4,916 ft, where the cattle with HAPH recovered, and then 
the cattle were bred within each phenotypic class. Bulls and 
cows were equally represented. First-generation offspring 
were studied at both low and high altitude. Mean PA 
pressures were 27±4 (STD] at low altitude. After residing 
at 10,000 ft for 2 months, the offspring of cattle with HAPH 
(brisket disease] had a mean PA pressure of 87 ±7 (SE], and 
those from resistant stock had a mean PA of 44±3 mmHg. 
Second-generation breeding (breeding experiments not 
overtly clear] within these two cohorts of calves yielded 
the same pattern of susceptibility and resistance. Pressures 
continued to diverge over time at altitude. Proof that the 
stimulus to pulmonary hypertension was hypoxia and not 
"altitude" was subsequently demonstrated by exposures 
in a hypoxic chamber at low altitude, whereby similar 
pulmonary hypertensive responses were recorded. Arterial 
blood gas measurements in susceptible and resistant cattle 
revealed similar severity of hypoxemia, similar alveolar 
ventilation (defined by the arterial carbon dioxide partial 
pressure, PC0 2 ] and only small differences in hematocrit, 
not considered sufficient to cause viscosity effects. The 
mode of inheritance of HPAH currently remains unclear. 181 
Nevertheless, what does remain clear among several studies 
is that HAPH in cattle is heritable, which indicates that cattle 
may be selected for resistance or susceptibility, and that the 
major gene(s] involved are likely to be discovered using 
large-scale genomic approaches followed by fine mapping. 

Our initial screen using four functional candidate genes 
(selected because they are involved in human PAH] was 
negative and suggested that a candidate gene approach 
may not be efficient. We then employed a 10K bovine 
SNP array for a GWA aimed at unearthing new potential 
candidate genes. This analysis revealed six or more 
genes of interest proximal to bovine SNPs exhibiting the 
most disparate distributions between severe HAPH and 
resistant cattle. Of the genes located near the SNPs of 
interest, NADH dehydrogenase (ubiquinone] flavoprotein 
2 (NDUFV], myosin heavy chain 15 (MYH15] and the 
myocardial signaling protein (FKBP1A] were candidates 
possibly involved in pulmonary hypertension. NADH 
dehydrogenase is no different between affected and 
resistant cattle, and additional studies are needed to 
further evaluate candidate genes. Furthermore, the recent 
availability of higher density bovine SNP arrays [44,45] 
provides a natural progression to more thorough GWA, 
which are likely to identify additional candidate genes due 
to greatly enhanced genomic coverage. The use of these 
new higher density bovine SNP arrays with larger cohorts 
of severe HAPH and resistant cattle has the potential to 
maximize the probability of detecting SNPs that display 
a nonrandom relationship with either severe HAPH or 



putative resistance. It is likely that the major gene(s] 
responsible for bovine HAPH will have identifiable human 
ortholog(s] that may also modify human pulmonary 
hypertension related to hypoxia (altitude, emphysema, 
hypoventilation syndromes], and possibly in other 
conditions that lead to pulmonary hypertension. 

Biological significance of HAPH signature 

To gain biological insights into resistance or susceptibility 
to HAPH, we analyzed the expression data using the Gene 
Ontology and Ingenuity Map database. Our goal was to 
identify the disease processes, physiological, cellular 
and molecular functions and the canonical pathways that 
are enriched in HAPH. These analyses showed that the 
top disease process reflected in our data is respiratory 
disease (Table 2]; the top three physiological processes 
are organism survival, hematological system and immune 
cell trafficking; the top three cellular and molecular 
pathways are cellular movement, cell signaling and small 
molecule biochemistry; and the top canonical pathways 
enriched in HAPH are IL-10 signaling, RXRactivation,and, 
interestingly, immune and endothelial cell function (data 
not shown but available on request]. 

We th en used the IPA application (IPA version 8.7] 
to examine the biological context and interactions of 
our signature. The IPA network analysis identified 10 
networks with high P-values (l.OxlO 21 ] (data not shown 
but available upon request]. The top three networks 
had many genes (MAP13, STAT5a/b, ERK1/2, Cyclins, 
VEGF, PDGF, IL1B, ID1, CREB, p38MAPK, IL15, etc.] 
with a documented role in endothelial and immune cell 
function through their effects on cellular pathways such 
as TGF-beta, BMP and MAPK (important in human PAH] 
and molecular functions such as cell death, survival, 
inflammatory response and cell morphology. Function 
annotation of all the networks showed that the top 
canonical functions associated with these 10 networks 
are BMP, TGF-beta, PPAR, IL6, endothelial cell function, 
NOS, TREM1, VEGF, glucocorticoid, ERK/MAPK, NF-kB 
and HIFla signaling (data not shown but available upon 
request]. Not surprisingly, many of these pathways 
are also known to be important in human pulmonary 
hypertension. 

In summary, herein, we have provided the results of 
the first molecular interrogation of HPAH, including a 
low-density SNP search, sequencing of four candidate 
genes and analysis of expression arrays to seek clusters 
of genes that may be a signature of hereditary bovine 
HAPH. Our hope is that this initial study of HPAH will 
spur further discovery and investigation of the molecular 
mechanisms behind HPAH. Our next challenge will be to 
attempt whole genome sequencing on the hypothesis that 
a polymorphism responsible for an autosomal-dominant 
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phenotype at altitude will emerge from the similar 
background genomes of the Angus herds. 
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